# autoplot(prcomp(t(mcp_zscore_matrix)), data = person_metadata_df, colour = 'person', main = "MCPCounter")
# autoplot(prcomp(t(mcp_zscore_matrix)), data = person_metadata_df, colour = 'cell_type', main = "MCPCounter")
# autoplot(prcomp(t(cs_zscore_matrix)), data = person_metadata_df, colour = 'person', main = "Cibersort")
# autoplot(prcomp(t(cs_zscore_matrix)), data = person_metadata_df, colour = 'cell_type', main = "Cibersort")
pheatmap(
    mcp_zscore_matrix, 
    main = "MCPCounter genes", 
    annotation_row = mcp_heatmap_row_df,
    annotation_col = heatmap_col_df,
    scale = "none")

pheatmap(
    cs_zscore_matrix,
    main = "Cibersort genes", 
    annotation_col = heatmap_col_df,
    scale = "none")

# pheatmap(
#     mcp_mean_zscore_matrix, 
#     main = "MCPCounter genes, collapsed by sample/celltype", 
#     annotation_row = mcp_mean_heatmap_row_df,
#     annotation_col = mean_heatmap_col_df,
#     scale = "none")
# pheatmap(
#     cs_mean_zscore_matrix,
#     main = "Cibersort genes, collapsed by sample/celltype", 
#     annotation_col = mean_heatmap_col_df,
#     scale = "none")
# kable(mcp_result_df)
# kable(cs_result_df)